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ABSTRACT 


Methods of investigating internal waves from bathythermograph 
information are discussed. A series of 651 hourly observations 
taken at Ocean Station "P" is analyzed by two methods to determine 
the characteristics of the internal waves present. One method 
involves a study of the fluctuations in depth of six selected 
isotherms; the other involves a study of the fluctuations in 
depth of the top, center, and bottom of the thermocline. Ап 
objective method is used to determine the depths of these 
features from the bathythermograph data. 

Spectra of these features are calculated by Fast Fourier 
transform, hanned, and smoothed, and these spectra are compared 
with previously published results obtained from a different 


method of analysis of the same time series. 
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I. INTRODUCTION 


Tidal period internal waves in the deep ocean were observed 
as long ago as 1907, by the Swedish oceanographer 
Pettersson 1909]. Little work was done in describing these 
waves or investigating their cause before 1950, although Stokes 
developed the theoretical treatment which could describe internal 
waves--neglecting geostrophic effects-- in 1847. This theory 
was used by Ekman [1904 in explaining the "dead water" 
phenomenon in Norwegian fjords. Although internal waves have 
been observed to exist for more than half a century, only 
recently have mathematical tools and equipment been developed 
which enable investigators to more closely analyze these waves in 
terms of periodic components and energy distribution. Although 
oceanographers are now in a position to say a great deal about 
the various frequencies which are present in a spectrum of 
internal waves, the presence of many of these frequencies 
and the causes of them have not yet been explained. 

The most useful tools for analysis of internal wave data 
are Fourier analysis and power spectrum analysis. Fourier 
analysis, although a tool of mathematicians for decades, 
required great expenditure of time to be carried out before the 
advent of high-speed digital computers, and has only recently 
been brought into use in analyzing internal wave records. Power 
spectrum analysis, developed and refined by Blackman and 
Tukey 1959] has been used by electronic and communications 


engineers for several years, but has only recently been applied 


to ocean wave analysis. Programming either of these methods for 
use on a digital computer makes possible analysis of large amounts 
of data accurately in a reasonably short time. 

In order to describe internal waves, it is at least necessary 
to obtain a time or space series of observations for analysis. 
The most common method of obtaining a time series is to observe 
temporal changes in thermal structure at a fixed location. This 
is sometimes done from weather observation ships located at deep 
ocean stations. One of these, Ocean Station "P", located in the 
northeastern Pacific Ocean at 145° West Longitude, 50° North 
Latitude and maintained by the Department of Transport of the 
Dominion of Canada, has been making bathythermograph observations 
at a maximum interval of 12 hours since July 1952. From time to 
time, bathythermograph observations are made from the ship 
occupying Station "P" at four hour intervals, or more frequently. 
The resulting bathythermograph data provides an excellent time 
series for analysis of internal waves. 

The data for use in this thesis consisted of hourly 
bathythermograph observations, obtained with mechanical bathy- 
thermographs from 2000 GMT, 5 August 1961 to 2200 GMT, 

l September 1961, a total of 651 hourly observations. The data 
were furnished in digital format on punched cards by the 
National Oceanographic Data Center, Washington, D. C. Digital 
analysis was performed on an IBM 360 computer at the Naval 


Postgraduate School Computer Facility, Monterey, California. 


Bathythermograms of the data used in this thesis may be found 
in Manuscript Report Series No. 106, published by the Fisheries 
Research Board of Canada, Programmed by the Canadian Committee on 


Oceanography. 


II. OBTAINING DATA FOR POWER SPECTRUM ANALYSIS 


One of the most commonly used methods of obtaining data for 
analysis in describing internal waves has been following the 
depth fluctuations of selected isotherms. This method was used 
by Tabata [1965] at Ocean Station "Р" іп Бе northeastern 
СТ" [1965 in the North Atlantic Ocean: 
Although this method has been quite successful in the past, its 
validity depends on rather specific conditions. Internal waves 
propagate along isopycnal surfaces, and have their greatest 
amplitude at a density discontinuity. An isothermal surface 
will closely approximate an isopycnal surface for propagation 
of an internal wave only if the water colum is nearly isohaline 
and the wave is of small amplitude. In this case one can approx- 
imate the change in density due to salinity and the change in 
density due to pressure fluctuations by zero. In choosing the 
isotherms of interest, Tabata selected isotherms which were 
within the permanent thermocline, assuming that the thermocline 
would oscillate as a unit, although phase аз лсе Е 
thermoclines would not be indicated in the power spectrum. 
Seiwell, in the North Atlantic, chose the 20°C isotherm as being 
close to defining the top of the thermocline, or depth of the 
mixed layer. 

Defant [1932] in examining the data collected by METEOR 
anchor stations, used a similar method, although he examined 
temperature fluctuations at fixed depths. In order to compare 


his results with that of other investigators, one must select 


a depth near the top of the thermocline, since other investigators 
used isotherms near the top of the thermocline for analysis. 
Decreases in temperature would indicate the approach of a crest, 
and increases in temperature the approach of a trough. 

LaFond and Rao [1954] obtain their raw data from the mixed 
layer depth, disregarding any changes above or below this depth, 
which corresponds to the top of the permanent thermocline. This 
would appear to be a more realistic approach, since the top of 
the thermocline will approximate a density discontinuity, and 
internal waves should have their greatest amplitudes at this 
point, and thus be most easily detected. 

In this thesis, the fluctuations in depth of the top, center, 
and bottom of the thermocline were used to investigate internal 
waves. The depths of these features were obtained using a method 
originated by Boston 1966. in which the thermocline is made to 
approximate a Gaussian distribution of depths and related 
temperatures. The depths thus obtained were subjected to a 
Fourier analysis, and the coefficients obtained were hanned 
and smoothed over discrete frequency bands. 

For comparison of results of isotherm analysis, the depths 
of six selected isotherms were followed through the time series 
and subjected to analytical procedures identical to those used 


in the spectral analysis of the thermocline. 


III. OBTAINING THERMOCLINE DEPTHS 


A. THE GAUSSIAN THERMOCLINE 


A gaussian, or normal, distribution of temperature (T) as a 


function of elevation (Z) is given by the equation 


Γη 


Т0) e “dz 


where Z is positive upwards. The frequency function corresponding 
to this distribution function is 


dTZ) 1 τἳ 


GZ ΠΠ 5 
Using a basic statistical approach, defining a normalized variable 


І DAR 
pa лс ги 


where m is the mean of the Z values, and O is the standard 


deviation of the array of Z values, we obtain a new frequency 


function which is 


GI. CN 2 
-а2 о/т 


This function has the following characteristics: 
1. It is symmetric about the point Z = m 
2. It has two symmetric points of inflection at Z = тї 
3. It has a maximum rate of change at Z = пс үЗ 
A change in the numerical value of m causes a displacement of the 


Curve in the vertical direction, but does not alter its form. 
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A change in ` has the effect of altering the scale in the 


+ 


vertical direction. The smaller the value of J , the more 


concentrated the curve is about the point Z m. 
If the water surface is defined to be Z - O, and the depth of 
the water considered to be infinite, then the function defining 
the vertical temperature distribution is 
(Z-mY- 
psu 
PER 
dZ 
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The center of the curve T(Z') versus Z is found at the point 
д = -m. 

If this curve is considered as a bathythermograph trace, the 
following definitions may be made. The center of the distribution 
of T(Z') versus Z is the center of the thermocline, called Z,. 

The top and bottom of the thermocline will correspond to the 
points of maximum rate of change of slope of this distribution. 

If Z, is used to designate the depth of the top of the thermocline 
апа 2ь+ to designate the depth of the bottom of the thermocline, 


these definitions become 


ue = -т 
Zt = -mO 
Zb t^ -m-O 


These features are illustrated in Figure 1А. 
To obtain the necessary values to apply these definitions, 
consider the first and second moments of the (T,Z) curve about 


the origin, Z - O. The nth moment of the curve about the origin 


LL 


is, by definition, the mean value of the array о called 2", 
where n is a positive integer. With a normal curve, the first 
moment is identical with the mean, which corresponds to the 


center of the thermocline. 


ТО) 


Zc: SICH ZT )aT 
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where T(O) is the temperature at the surface and T(-9O) is the 
temperature of the bottom of the lower isothermal layer which, 
so far, has been assumed to be the bottom of the ocean. The 


second moment 


T(O) 
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τος 


can be used to compute O , since the variance of the distribution, 


τ 


J i , is defined as 
ots 87 

In practice, several approximations to this exact solution 
must be made. The temperature T(-.) is measured at the bottom 
of the bathythermograph trace. The differentials, dT, become 
finite dup cq AD and the integrals are replaced by 
summations. 

In changing from an integral representation to a summation 
over finite intervals, a certain amount of error in the value 
of ." is introduced, due to the fact that summation in this 


fashion is not a continuous process, whereas integration is. 


JE 


To compensate for this loss of accuracy, a weighting factor must 
be introduced in all computations involving О’. Grosfils [1968 


calls this factor sk, and bases the size of the factor on the 


EN 
saze or OT. 


For OT21.0° sk = 1.30 
For Ӧт= 0.59 sk = 1.47 

Intermediate values of OT will have intermediate values of sk, 

but these values were chosen as being those most commonly 

encountered. Using this factor, the equations defining the top 

and bottom of the thermocline become 


rt = -m + O sk 


bi -m - O sk 


Using this approach to define the thermocline, one can fully 
describe a thermocline knowing the temperature at the surface, 
the temperature at the bottom of the bathythermograph cast, and 
depths corresponding to some constant temperature interval. 


The center of the Gaussian thermocline can now be defined as 


ia 1 | 
Z-= - Ini Т; э ) АТ.) 


апа ѕіпсе 
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then 


Using these factors, the top and the bottom of the Gaussian 


thermocline may readily be determined. 
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In the event the distribution of depth with temperature is 
not Gaussian, the inaccuracies built into this objective 
definition are intolerable. Whether or not this distribution 
is Gaussian can be determined by computing the third moment of 
the (T,Z) curve about the origin. This, however, is unnecessary 


when one applies a non-Gaussian approach to the thermal structure. 


B. THE NON-GAUSSIAN THERMOCLINE 

The approach that was used by Boston and by Grosfils was 
based on depth data given at constant temperature intervals. The 
data used for this thesis were temperature data given at 
constant depth intervals. This difference does not affect 
the treatment of the data except for the determination of the 
temporary center of the thermocline. The details of the method 
used in this thesis to find the center of the thermocline are 
given in Section IV. 

When the depth of the center of the thermocline is found, by 
whatever method the data requires, the data is then processed as 
two Gaussian thermoclines--one based on the depth distribution of 
temperatures above and including the depth Z,, the other based оп 
the depth distribution of temperatures below and including the 
depth Zg. 

Іп order to locate Z,., onily the upper portion of the (T,Z) 
Curve is considered. This portion of the curve is "reflected" 
about Zo, to produce the curve shown in Figure 1B. From this 
curve (-- is Computed by methods previously described and, 


using this value, Z+ may be calculated. 
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In order to determine Zpr, a similar procedure is followed, 
except that in this case, the lower portion of the (T,Z) curve 
is "reflected" about Zœ, to produce the curve shown in Figure lC. 
Note that it is possible for this resulting curve to include 
negative depths, but since the depths are squared in 
computing сло this does not result in a narrowing of the 
distribution through an underestimate of O . 

It can be seen that in the event the thermocline is Gaussian, 
treating it in this fashion will not alter the results obtained 
by treating it as Gaussian. Thus, every(T,Z) input may be treated 
as a non-Gaussian distribution, regardless of its true form. 
This is the method which was used in processing the bathythermogram 


data for this thesis. 
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IV. DIGITAL METHODS USED TO OBTAIN DEPTH INPUTS 


The bathythermograph data was obtained from the National 
Oceanographic Data Center in digital format. For each 
bathythermograph trace, the digital data was contained on three 
or four computer data cards, including the sequential number of 
the observation, the date and hour of the observation, 
temperature in degrees Centigrade at each five meters of depth, 
and such other bits of synoptic information as the location of 
the observation, the ship making the observation, and a number 
identifying the cruise during which the data was taken (Figure 2). 

The first phase of data reduction was to transfer the 
sequential number of the observation, the time and date, and 
the temperatures at the first twenty depths--which correspond to 
temperatures from the surface to a depth of 95 meters--to a 
computer data card. This was not essential for data reduction, 
but was done for convenience as it reduced the size of the data 
deck by three quarters. By inspection of bathythermograms for 
the entire series of observations the permanent thermocline was 
in all cases determined to be shallower than 95 meters in depth. 

The data deck thus obtained was then analyzed to obtain the 
depths of the top, center, and bottom of the thermocline using 
a method based on a digital computer program developed by 
Grosfils |1968 . Grosfils' program, however, was written for 
data consisting of depths corresponding to equally spaced 
temperatures, and since the NODC format differs from this, some 


modifications had to be made in the program. 
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The computer first reads the data from the data deck, then 
calls subroutine INTERP, which fits straight lines through the 
given data points, then determines the depths of temperatures at 
intervals of 0.1°C from the surface to 95 meters. The temperature 
is read by the computer as an integer--in tenths of degrees--and 
left in this form for the interpolation process, to eliminate 
truncation error in manipulating the temperatures. 

In the program written by Grosfils, the center of the 
thermocline was defined as the point where the second difference 
of depths was greater than Zero. Due to the format of the data 
used in his program, however, a great deal more information was 
known about the thermocline than is given in the NODC format. 
For this reason, the following method was used to define the 
center of the thermocline, utilizing the array of temperatures 
before interpolation. 

In subroutine INTERP, the first differences of temperatures 
were calculated, from the surface (corresponding to T4) downward. 


The ith first difference is defined as 


The first time the ith first difference is less than 29, 3C ; т. 15 
stored as ITTOP and the computer then begins calculating first 
differences again at a depth of 45 meters. When the ith first 
difference then becomes greater than ο σα the computer stores 
T;,+j as ITBOT, The temperature of the center of the thermocline 
is then defined as the arithmetic mean of ITTOP and ITBOT, and is 
stored as ITC. The computer then interpolates depths and 


temperatures as descrived previously. 
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Once the interpolation is completed, subroutine PURGE is 
called to insure that only one depth is recorded for each 
temperature value, which would not be the case if any temperature 
appeared more than once in the data which was read in. PURGE 
examines the interpolated temperatures, and in the event that 
any temperature is greater than or equal to a temperature at a 
shallower depth, discards the corresponding depth. After 
proceeding through the entire interpolated trace in this manner, 
the trace is reduced to one in which there is only one depth 
corresponding to any given temperature, and it will be the 
shallowest depth at which this temperature appears. Although 
the program at this point has destroyed all information about 
inversions, this was not felt to be an unreasonable step. The 
raw data were all examined for inversions prior to the interpo- 
lation step, and in no case were inversions present greater than 
40.1 C in five meters of depth. 

At this point, subroutine FINDZC is called, and it examines 
the temperature and depth information thus computed and chooses 
the depth of the first (shallowest) temperature which is less 
than or equal to ITC. It stores this depth as ZC, the depth of 
the center of the thermocline. This depth, along with all the 
remaining temperature and depth information, is returned to the 
main program. 


Using this information, the main program calls subroutine 


= а 9. б 
у 1 H 3 ` CN 1 % 3 


NONGA, which uses the procedure described in Section II to find 
the top and bottom of the thermocline through the use of 


subroutine GAUSSI and GAUSS2, respectively. 
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For comparison with Tabata'!'s results, the depth fluctuations 
of selected isotherms were also subjected to Fourier analysis. 
Depth information for this analysis was obtained in a manner 
similar to that described above, except that for depths of 
isotherms, subroutine FINDZS is called from INTERP after PURGE 
and the depths are determined in a manner similar to that used to 
find ZC, 

Fluctuations of isotherms and thermoclines with time are 


shown in Figures 3 and 4, respectively. 
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V. COMPUTATION OF POWER SPECTRA 


The reworked data obtained as described in Section IV were 
subjected to Fourier analysis to produce a power spectrum. 
Subroutine RHARM, which is contained in the IBM Scientific 
Subroutine Library, was used for this analysis. This subroutine 
uses an algorithm, developed by Tukey and Cooley 1965), to 


transform the depth fluctuations to a Fourier series of the form 


7 = \ COST Nares Æt 

E 

n=O 
where Wis the frequency of a given sinusoidal component in 
radians per hour, N is half the total number of data points plus 
one, t is time in hours, and /] is the instantaneous wave elevation 
above its mean position. The sum of the squares of the 
coefficients, асты е^, is proportional to the power contained in 
the nth component, whose period is N/(n) hours. ag is the mean 
value of the data, and to eliminate leakage of this large value 
into neighboring frequency components, all data was adjusted to 
a mean value of zero in subroutine AVERGE before RHARM was called. 
bg is always equal to zero, and is therefore ignored in subsequent 
manipulations. 

Once the Fourier coefficients were obtained, they were hanned 
by subroutine HANN. This subroutine used the following recursion 
formula to refine the coefficients: 

- 0.5а, - 0.25а 


1 2 
а. =-0 „25а а + 0.5a; = 0.25a;+1 i= 2,3,...,N-1 


d 


+ fe о р > 


ак 7 0-25ац-1 + 0.5ак 


The sine coefficients were similarly treated. Hanning the data is 
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a method of greatly reducing leakage of power from frequencies 
containing large amounts of power into neighboring frequencies. 
According to Bingham, et. al. |1967 | if the Fourier transform is 
perfect forw,, Wo, W3,... and W is none of the Wi, the 
leakage without hanning is proportional to |а |2 and hanning 
the coefficients reduces the leakage to a term proportional 

to hær. 

The hanned data was then plotted on a graph of power (spectral 
density) versus frequency. Although this presentation gave 
indications of the frequencies of the major components of the 
internal waves present, it was objectionably noisy. To reduce 
this noise, a method of power spectrum smoothing was utilized to 
produce the final spectrum. This method, originated by 
Jones 1965 was used in subroutine SMOOTH. The method consists 


of first arranging the squares of the non-zero a and b coef- 


ficients in order and renaming them s, as follows: 


2 2 2 2 2 
ag »8ι ‚Бу ‚а> T sS. 


5] 52 534 25η э56 »..... 
The smoothed spectrum is obtained by averaging these numbers in 
groups of 2m, except for the first group, which contains m, and 
the last group, which contains those numbers which do not fit 
into groups of (2p * l)m, where p is the number of groups of 


em included in the data. 
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Зо = N/m кл T e Р 
3 7-76; 
51 = N/m Si 
me 1 
- m 
(2р1 


where mp is the number of pairs remaining in the last step and 
N is the total number of Fourier coefficients used in the 
smoothing. By this method of smoothing, Sọ will be centered at 
.25/N ^t, 5] will be centered at YN Ot, 52 аї 2/қ Зе, and so 
forth. Care must be taken to choose m sufficiently large that 
significant smoothing takes place without choosing it so large 
that it eliminates detail from the final spectrum. 

Subroutine SMOOTH was run with m's of 2, 3, 4, 5, 10, 20, 
and 50. Of these, the best results were obtained with m= 4, 
and this grouping factor was used in the final analysis. After 
the data were smoothed, subroutine PERIOD was called to generate 
counters to be used for the final graphic output. 

The data were broken into two groups, the first containing 
the first 512 data points (consecutive numbers 1 through 512) and 
the second containing the last 512 data points (consecutive 
numbers 139 through 651). More confidence may be placed in 
spectral estimates if the spectral estimates are the averages 
of several sets of independent data. The nearest approximation 


to this which was possible with this data was the treatment which 
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was used. Had a sample size of 256 points been chosen, totally 
independent data sets could have been used, and confidence in 

the results would have been increased. The reduction in the 
number of data points by a factor of two, however, would have so 
significantly reduced the accuracy of the harmonic analysis that 
the increased confidence in the results would have been 
meaningless. For this reason, the use of two partially independ- 
ent data sets of 512 points each was decided upon. 

The two groups of hanned, smoothed data were then added 
together, term for term, by subroutine GROUP, and averaged. 
Plots of the results of this averaging for the top, center, and 
bottom of the thermocline are shown in Figures 5, 6, and 7. 

Analysis of the depth fluctuations of selected isotherms was 
carried on in the same manner as analysis of the thermocline 
data. The smoothed spectra of fluctuations of the isotherms 
analyzed are shown in Figures 8, 9, 10, ll, 12, and 13. Рог 
comparison, Tabata's results for the 6°C and 10°C isotherms are 
shown in Figure lH. 

The complete computer program used in data reduction is 


given in the Computer Program section of this thesis. 
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VI. RESULTS 


In his paper, Tabata reports a strong semidiurnal component 
in the internal waves analyzed through depth fluctuations of the 
6°C through 10°С isotherms, with greater power in higher 
frequencies noted in analysis of the deeper isotherms. This 
general behavior was also noted in the digitally analyzed data 
used in this thesis. 

A power peak centered at an 11.64 hour period was present 
in the spectra for all the isotherms analyzed, as well as in the 
spectra of the top, center, and bottom of the thermocline--the 
shallowest feature studied--also shows a peak corresponding to a 
period of 25.6 hours, probably associated with a principal lunar 
diurnal (0)) tidal component. 

The power spectrum due to fluctuations of the bottom of the 
thermocline indicates a shift of power to higher frequencies, 
with peaks associated with periods of approximately 4.5 and 3.5 
hours, and a maximum energy associated with waves of 9.14 hour 
period. Tabata reports waves of periods equal to five hours 
and slightly less than four hours present for the deepest of 
the isotherms which he studied, with which these data agree. In 
no case, however, does he report waves of periods near nine 
hours. Relative peaks in spectral density having periods between 
nine and ten hours (most often at a period of 9.85 hours) were 
found associated with all features. Waves with periods of 
16.00 hours, which is very close to the inértial :peridd iof 15.7. 
hours at Ocean Station "P" were also noted in all spectra 


produced. 
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In the original plan for the analysis, the fluctuations of the 
thickness of the thermocline were also to be investigated, to 
determine whether the thermocline was oscillating in the funda- 
mental mode--in which the entire thermocline oscillates as a 
unit--or mode two--in which the top and bottom of the thermocline 
are 180° out of phase, resulting in a periodic thickening and 
thinning of the thermocline's vertical extent. The spectra pro- 
duced by analysis of the thickness showed nearly constant spec- 
tral density throughout the entire range of frequencies studied, 
except for large contributions in the low frequency end of the 
Spectrum. This suggested that the thermocline was oscillating in 
the fundamental mode, and the spectral densities at high 
frequencies were due to noise. To verify this, the correlations 
of the top, center, and bottom of the thermocline were computed. 
It was found that the depth fluctuations of each of these 
features with the others was greater than +.90 which confirms 
the fact that the entire thermocline was oscillating as a unit. 
This conclusion agrees with the results published by Tabata, who 
found phase differences on the order of i 20? between oscillations 
» of the isotherms he analyzed. 

Ihe spectrum associated with the top of the thermocline 
Showed a very large contribution due to waves of very long 
period. Examination of Figure 3 shows the reason for this. The 
top of the thermocline was slowly sinking with time, presumably 
due to mixing in the upper layer. This effect, although small 
in magnitude, would have been identified in the analysis as an 


internal wave with a period equal to the total length of the 
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record, which in this case was 512 hours. In addition, during 
the period 16 August to 21 August, the water column, which was 
almost perfectly isothermal above and below the thermocline at 
all other times, developed near-surface transient thermoclines. 
As the non-Gaussian analysis used to obtain thermocline depths 
did not eliminate transient thermoclines before computing the 
depths associated with the main thermocline, these near-surface 
transients resulted in the calculation of an abnormally large 
variance by subroutine GAUSS1, which had the effect of moving 
the calculated top of the thermocline closer to the surface than 
it actually was. The fact that the time period during which this 
occurred was included in both sets of 512 hourly observations 
meant that an additional large power component associated with 

a period of 512 hours was found in the analysis. This is an 
apparent effect, rather than a real one, and even when combined 
with the long-term sinking of the thermocline, fails to obscure 
the results in the frequency range of interest. 

Examination of Figure 3 also reveals that the center of the 
thermocline also showed a slight downward displacement with time, 
but of less magnitude than that of the top of the thermocline. 
This resulted in a relatively large amount of power associated 
with long-period internal waves at the center of the thermo- 
cline, but much smaller than the power contribution to the 
spectrum of the top of the thermocline. The bottom of the 
thermocline showed no change in mean depth with time, which 
indicates that surface layer mixing extended only to a depth 


Slightly below the center of the thermocline. Although in most 


cases large amounts of power were associated with low frequencies, 
these frequencies were not considered of interest, and are not 
indicated on the graphs of the final spectra. A complete presen- 
tation of all dra produced by the analysis is presented in the 
COMPUTER OUTPUT section. 

Tabata's results for the 6?C and 10?C isotherms are presented 
in Figure lH in a semilog plot for comparison with the results oí 
the digital analysis. Although these curves differ in detail 
from those presented for the same isotherms analyzed digitally in 
this thesis, the overall form is strikingly similar. For both 
analyses of the 10°C isotherm, the greatest peak spectral density 
is associated with waves whose period is greater than the semi- 
diurnal tidal period. For both analyses of the 6°C isotherm, 
the spectral density is rather well-distributed throughout the 
higher frequency components, which also have a greater mean spec- 
tral density than the higher frequency components of the analyzed 
10°C isotherm data. This indicates that the higher frequency 
components showed an increase in amplitude with increasing depth, 
a behavior which was also noted by Haurwitz, Stommel, and Munk 
1959]. 

The Brunt-Vaisala frequency associated with the thermal 
structure of this time series is 3.79 x 1072 зес"1. This indi- 
cates that the thermal structure is capable of supporting inter- 
nal waves of any period greater than 26.4 seconds. The Nyquist 
frequency for a digital time series sampled at hourly intervals 
is 0.5 Һет, Thus, there is a large range of frequencies which 


may be present but of which the data does not permit investigation. 
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The causes of the internal waves with periods near 24, 16, 
or 12 hours seem obvious. The 24 and 12 hour periods are 
probably associated with internal waves which are being driven 
by diurnal or semidiurnal tidal forces. As mentioned previously, 
the 16 hour period seems to be associated with inertial iinternal 
waves. The causes for the internal waves with periods between 
four and six hours--indicated both in this thesis and in Tabata's 
study--are not so easily theorized. 

Another possible source of internal waves is baroclinic 
Rossby waves, but this type of wave is associated with extremely 
low frequencies, with periods on the order of days, so it would 
seem that these waves will not explain the higher frequency 
waves indicated by the analysis. Two other possible causes for 
these waves do exist, however. They might be long-period waves 
which have been reflected from a coast and have interacted with 
varying density gradients, which act as a weak filter, but this 
theory cannot be tested unless one knows the direction of 
propagation of the higher frequency internal waves. 

The other possible source of these waves is a non-linear 
interaction between surface pressure.disturbances, but in order 
to examine this as a possible cause, one must have records of 
both surface waves and weather. Although weather records exist 
for the time during which the bathythermographs were taken, no 
surface wave information is available. To be of value in 
investigating this process as a causative agent, the surface 
wave information would have to be very detailed, and wave records 


of the required detail can only be obtained in the deep ocean 
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through highly sophisticated and careful measurements which are 
not compatible with the mission of a weather ship. 

In conclusion, it seems that little can be said about these 
higher-frequency internal waves except that they exist. In order 
to say more, one would require a great deal more data than is 
available. Little has been done in the investigation of internal 
waves to date, and one can see that it represents a fertile 


subject for further research. 


43 


VII. RECOMMENDATIONS 


During the execution of this thesis, a number of changes in 
the format and amount of data were wished. Had these changes 
been made, the analysis of the data would have been simplified 
and the accuracy of the results enhanced. For these reasons, the 
Following recommendations are made. 

In digitizing mechanical bathythermograph traces, it would 
greatly assist the user of the digitized data if the data 
through a thermocline were recorded as depths at equally spaced 
temperatures, instead of temperatures at equaily spaced depths 
as is now done. This would provide an increased number of points 
in the region of maximum gradient, enabling the digitized data 
to be used to more accurately reconstruct the original trace. 

The method presently used to digitize expendable bathythermograph 
traces, which records a data point at the location where the 
trace departs from linearity by a specified amount, allows use 

of digitized data to almost perfectly reconstruct the trace. 

If salinity data were also available for analysis, it would 
be a simple task to compute depths of isopycnals using a 
computer. If terms to second order involving salinity and 
pressure were introduced to the already-known temperature 
information, depth-density profiles of much greater accuracy 
would result, and this in turn could be analyzed to produce 
internal wave information of higher quality. 

Closer spacing of observations would enable analysis of 
higher frequency components, since through aliasing the smallest 


wave period that can be detected using a sampling interval of one 


цц 


hour is two hours. Also, as was suggested by Tabata, utilizin, 
data presently available over extremely long periods, with t as 
large as a year, would allow analysis of extremely low frequency 
internal waves. 

As was indicated in the Results section of this thesis, 
there is a need for information regarding the direction of 
propagation of the internal waves indicated. This could be got- 
ten using a cross array, consisting of two orthogonal arrays 
ot the type described by Gilchrist 1966) For use in investi- 
guting the direction of propagation of surface waves. Several 
such arrays would be required to allow investigation of a large 
range of possible wavelengths, but it would seem that the 
information obtained as a result of such an expenditure of effort 
and equipment would be a highly significant contribution to the 


present knowledge of internal waves. 
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APPENDIX A 


EXAMPLE OF COMPUTATION 


The following four figures indicate graphically the operation 
of the subroutines used in computing the depths used in harmonica 
analysis. The data is from bathythermograph number 392 (Fig. 2).: 

Figure 15 indicates the operation of subroutine INTERP in fit- 
ting straight line segments between raw data points, and the 
choosing of a temperature corresponding to the center of ther- 
mocline (ITC). 

Figure 16 shows the operation of subroutine PURGE in elimi- 
nating small inversions and duplications of temperatures, reducing 
the data to a set of temperature-depth pairs from the surface to 
the bottom of the trace at an interval of 0.1°C. This figure 
also indicates the operation of subroutines FINDZC and FINDZS in 
determining the depth of the center of the non-Gaussian thermo- 
cline and the depths of the selected isotherms. 

Figures 17 and 18 show the thermal structures used in GAUSSI 
and GAUSS? to obtain the depth of the top and bottom of the 


thermocline, respectively. 
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SPECTRAL ESTIMATES FOR GROUP NUMBER 


19 DEG S DEG 


11 СЕС 


PERIOC 


М СО СЉ ОС (NJ O vf .O vf O. tr* O — WY G try C) rry f U vf <O CNJ VA OO TO DJT =T ~ OMN OM OOD LK OD Sr OA LHR чо ч ОО оС Сч З 
че ДСО АГ ry = QA T. UA CNJ GO CO UTA Nf SOU OL у =ч σώσω οσο CO Р UN WA O BHO FAINT NOW CIM c IPS xT O CO NY COUP Sun 
( Ө о @ Ө Ө о @ 9 9 € Ө 9 9 09 Ө 9 9 9 09 06 09 9 9 9 9 0 9 0 9 9 0 Ө 9 9 60 Φ 0 9 Ө 9 9 9 9 9 0 06 0 Ө 6 @  @ ° P Ө Ө 9 ә 6 9 0@ 6 9 9 9 ө 
СОШ ОЗ МИ +f vf CO CO (NJ O CO UJ O tü rN. vf (N GN] CO = ST CO FT CO = 00 vf LY CO ~= wN ANT O TiN О О чо ЧО ОО ОУ чо рч орт ОС У се 
ОМ у ЧСИ ч ΩΩ ΟΩώΨ(ΩΟΩΟΟΟΟ ¿ry CO ON MT, eo emm mex ны - ο ον ч ον κ 

= 


O rn CO QO. r— r= uy (N O Q — cr r= — cN un .O UNAY uy rn (N QO СЫ —r— ООУ С СООЗ ОСУ ООЗ ОФ SA —4 7 f OO NT e PSU CO GQ = ОМ ч «γω 
ОМ СОТ ООММ Ue OO C CQ (3 O O O CG: (N — ч ОсОО С О О СУ QO UA (NO CO +f (N O OQ e xr m = (N LIYA O OO (N O —— + — CO + (N C uA УЗ 

e o ө ө ө ө ө ө ө ө 9 9 ө ө ө ө ө ө ө ө ө 0 9 0 ө ө ө ө ө 0 ө ө ө 99 Φ 9 ө 99 9 9 о ооо оо о оо оо о оо о о ө ° ө ө ө ө 
и ры SE смао о us 
SQ 0 49 00 4 AT CN tr) AMAN OTOMONAN προ σπορ я = a xt a 


NO OD (0 00 ©) OD WI DD CM OES OO ex U^ C9 С rry J U" WA GJ GO =4 CGO QD uU P9 MDM MIND SIO CO (NU) = 4T. e 0090 O CN O^ uU 0 CO NQ U^ ND em m US C0 ait © 
=O TAL MEIN OO DIN OM xT VÀ = (N CO r= CN O Q^ 00 e —4 UN sf O UN tO CO (NJ CN = $r CO —4 CQ `Q ООО СОМ ч со ОО Y m QO O СО T (N= — — “O = 
6 6 о 0 оо о о о о о о о 6 0ὁ 90 Ө 90 9 9 0 9 90 9 0909 9 09 о о о 9$ @ о Φ 6 @ 6 Ө 9 9 9 09 09 0609 9 9 9 е 6 Ө Ө 9 0 о 9 ев 9 ө ө 
юр ыш ыы s 
BINAN = FNAT wO «АСОМ COM, Cem ON em mi «ἡ =i BE і 


C) CD OQ r ССГ C) CN O S nm ur O09 00 exa OQOQCON t= се (N CN Sl MOOV OV OPlP CO C2 MA CO e xf CO CN P9 mL ND e ОСУГ fr US un e P rry CO < rn Un" 
CIM OOO OM NON DOO DW BAO DN S exO0uYn mg rud oe O GO r= Q ú Jf (Q (N CN = Q GQ" О OSPF О ОЧАИ T V ANANN CN a e і С) СС 
о о Ө Ө Ө @ @ ϐ 0 Ө 09 0 606 Φ 9 Φ Φ Φ Φ Φ Φ Ө 9 9 9 9 9 09 0 9 9 9 9 Ө 9 0 9 9 9 @ 9 0 0 9$ е 0 9 0 9 0 6 9 0 9 9 9 9 99 Ө Ө 0 9? 9 @ 
|NO со + CNONUS e 20 4)O XC e G; C" G: O CQ i= = sO Q QO VA UA @JN UA $ + +f +f f $ y CY rO (ry rn (n cnn eo en Соу сеу СЧ СМ СУ СЧ СМ NN NN NAN NNN NN NNN NN ч 
«СМ MT (O (NJ UN =—4 =4 w чі 4 

CJ e і 


55 
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